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ABSTRACT 

We report first results from an ongoing monitoring campaign to measure time delays between the 
six images of the quasar SDSS J2222+2745, gravitationally lensed by a galaxy cluster. The time delay 
between A and B, the two most highly magnified images, is measured to be tab = 47.7 ± 6.0 days 
(95% confidence interval), consistent with previous model predictions for this lens system. The strong 
intrinsic variability of the quasar also allows us to derive a time delay value of tca = 722 ± 24 days 
between image C and A, in spite of modest overlap between their light curves in the current data 
set. Image C, which is predicted to lead all the other lensed quasar images, has undergone a sharp, 
monotonic flux increase of 60-75% during 2014. A corresponding brightening is firmly predicted to 
occur in images A and B during 2016. The amplitude of this rise indicates that time delays involving 
all six known images in this system, including those of the demagnified central images D-F, will be 
obtainable from further ground-based monitoring of this system during the next few years. 

Subject headinqs: galaxies: clusters: general — gravitational leasing — quasars: individual 
(SDSSJ2222-h2745) 


1. INTRODUCTION 

Quasars gravitationally lensed into multiple images by 
foreground clusters of galaxies are exceptionally rare. 
The two first reported cases were SDSS J1004-I-4112 (In- 
ada et al. 2003) and SDSS J1029-I-2623 (Inada et al. 
2006). A third case, SDSS J2222-I-2745, was discovered 
by Dahle et al. (2013), displaying six images (of which 
five were spectroscopically confirmed) of a quasar at red- 
shift Za = 2.82, with a maximum separation of 1571, 
lensed by a foreground cluster at Zs = 0.49. In con¬ 
trast, the Master Lens Database of known gravitational 
lens systems (Moustakas et al. 2012) lists ^ 120 known 
galaxy-scale quasar lenses. 

Refsdal (1964) showed how the relative time delays 
measured in multiply lensed quasar images can be used to 
measure the Hubble parameter Hq. In addition, time de¬ 
lay measurements can probe the underlying cosmological 
model (Refsdal 1966), yielding dark energy constraints 
which are complementary to other such probes (Linder 
2011; Treu et al. 2013; Sereno & Paraficz 2014). Time 
delay measurements of individual systems are subject to 
a ’’cosmic variance” of several percent, caused by the 
mass fluctuations of intervening structure along the line 
of sight (e.g., Bar-Kana 1996; Wambsganss et al. 2005). 
Aside from the cosmographical applications of time de¬ 
lays, they also probe the surface density distribution of 
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the lens, over the range of radii at which images are 
seen and time delays are measured. A unique property 
of SDSS J2222-I-2745 is the availability, at small clus¬ 
tercentric radii, of three fainter images (D-F) which are 
predicted to be demagnified, according to the six-image 
model presented by Dahle et al. (2013). These three 
images are still sufficiently bright to have their time de¬ 
lays measured from ground-based data, given the level of 
quasar variability reported in this paper. 

The longest time delays produced by cluster-scale 
lenses can be of order several years (Fohlmeister et al. 
2008; 2013), rather than weeks or months in the case 
of galaxy-scale lenses. With typical monitoring cadence 
of a few days to a few weeks, and given the measured 
trend of quasar variability amplitude to increase as a 
function of time scale (at least up to time scales of several 
years; e.g., Vanden Berk et al. 2004), cluster-scale lenses 
may yield smaller fractional uncertainties on the time 
delays, making cosmic variance-limited time delay mea¬ 
surements easier to obtain than for galaxy-scale lenses. A 
further advantage of cluster-scale lenses is the availability 
of additional constraints on the lensing mass distribution 
from other multiply lensed background galaxies as exem¬ 
plified by the HST-based studies of SDSS J1004-I-4112 
(Sharon et al. 2005; Oguri 2010) and SDSS J1029-b2623 
(Oguri et al. 2013). 

Measuring the time delays in wide-separation cluster- 
lensed quasars also enables three-dimensional studies of 
outflowing winds within the quasar itself by probing mul¬ 
tiple sightlines, as shown by Misawa et al. (2014). In 
this case, the time delays need to be known in order to 
disentangle time variability of absorption profiles in the 
quasar image spectra from genuine differences between 
the sightlines. 

Finally, knowing the time delay is crucial for advance 
planning of monitoring with increased frequency in antic¬ 
ipation of predicted strong flux changes, e.g., to measure 
the time delays of the faintest images in the system or 
to conduct reverberation mapping studies to measure the 
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mass of the black hole that powers the quasar (Blandford 
& McKee 1982; Peterson 1993, 2014). 

In this paper, we summarise results obtained from the 
first three years of an ongoing monitoring campaign of 
SDSS J2222+2745. In §[2]we present the data set avail¬ 
able so far, which is used to derive a value of the time 
delay tab in § [3] and a value of the longer delay tqa in 
m In§El we summarise and discuss the results, includ¬ 
ing a strong rise in the flux of image C during 2014, and 
the resulting prospects for future studies of this system. 

2. DATA 

A monitoring campaign to measure time delays in 
SDSS J2222-I-2745 was initiated at the 2.56m Nordic Op¬ 
tical Telescope (NOT) in October, 2012. The results 
in this paper are based on data from this ongoing cam¬ 
paign, with the most recent data obtained in August, 
2015. The NOT monitoring program makes use of the 
Andalucia Faint Object Spectrograph and Camera (AL- 
FOSC) instrument, which has a 2048^ pixel CCD detec¬ 
tor with 32-bit intensity resolution and a pixel scale of 
07188 pixel'At each monitoring epoch, a set of either 
3 X 300s or 3 X 600s exposures were obtained. 

The annual visibility period of SDSS J2222-I-2745 ex¬ 
tends from late April until early January. We chose the 
SDSS 5 -band for our monitoring campaign, in order to 
maximise the signal-to-noise of the quasar images, rel¬ 
ative to the galaxies in the 2 ; = 0.49 lensing cluster on 
which the central D-F images are superposed. This band 
is well centered on the prominent Lyman-alpha line in the 
spectrum of the 0 = 2.82 quasar, at an observed wave¬ 
length of 4640 A. 

In addition to the ALFOSC data, we include 5 -band 
data obtained as part of the original discovery of this 
lens (for further details, see Dahle et al. 2013), using the 
MOSaic CAmera (MOSCA) at NOT. MOSCA is a 2 x 2 
mosaic of four 2048 x 2048 CCDs used in 2 x 2 binned 
mode, yielding a pixel scale of 07217 pixel”MOSCA 
5 -band imaging was obtained for three separate epochs, 
on 2011 Sep. 24.93, 2012 Sep. 12.96, and 2012 Sep. 15.98 
UT, respectively. We also include the original SDSS DR 8 
5 -band photometry for this system, obtained on 2009 Oc¬ 
tober 18.19 UT. Finally, this paper also includes recently 
acquired Gemini imaging of SDSS J2222-I-2745 (taken 
under program GN-2015A-FT-16); these few images rep¬ 
resent the beginning of a new observational campaign to 
measure time delays for the much fainter central D-F 
images. Those data were acquired in 5 -band with Gem¬ 
ini North, using the GMOS instrument in imaging mode, 
on 2015 May 25.54, 2015 May 26.52, and 2015 July 25.48 
UT. 

In this paper, we report photometric measurements of 
the three brightest (A-G) images of SDSS J2222-I-2745. 
Unlike the fainter (D-F) images, these are all well sep¬ 
arated from foreground objects and their fluxes can be 
reliably measured without modelling and removal of ad¬ 
ditional sources. In order to calibrate the observed fluxes, 
we use a set of 5 reference stars (S1-S5; see Fig. [T|). In 
order to minimise the effect of spatial variations of the 
PSF, these are chosen within a radius V from the quasar 
images. In addition, the reference stars have blue col¬ 
ors similar to that of the quasar and do not have any 
nearby objects which may affect the flux estimate. Even 
during the best seeing (FWHM ^ 075), these stars have 


measured peak fluxes < 80kADU, well within the lin¬ 
ear regime of the ALFOSC detector, which extends up 
to ~ 350 — 400kADU. The photometric measurements 
in our data were all based on aperture photometry of 
the quasar images and reference stars, using the same 
aperture. An aperture correction was implicitly made by 
setting the magnitude zero point of each exposure such 
that the mean magnitude value of the ensemble of 5 ref¬ 
erence stars equals their mean 5 -band magnitude value 
in the SDSS (5 = 19.092). The accuracy of the time 
delay measurements reported in this paper is not depen¬ 
dent on a highly accurate absolute flux calibration, as 
long as the reference stars do not vary significantly in 
brightness. The validity of using S1-S5 as flux reference 
sources was tested by measuring the rms fluctuation in 
the magnitude offset of each star from the ensemble aver¬ 
age, over the time span covered by our ALFOSC observa¬ 
tions. The rms value of these fluctuations range between 
0.0061 and 0.0095 magnitude. This is two orders of mag¬ 
nitude smaller than the amplitude of the flux variations 
measured in SDSS J2222-I-2745, as described below. 

The ALFOSC imaging data were de-biased and flat 
fielded using standard techniques. Flat-fielding was typ¬ 
ically performed using the median of a series of flat field 
frames of the twilight sky obtained during the same night 
as the photometric data. As the ALFOSC detector is cos¬ 
metically clean, the three individual exposures for each 
monitoring epoch were treated individually in our anal¬ 
ysis rather than combined into a single frame. 

The Gemini GMOS imaging data were de-biased and 
flat fielded, again using standard techniques, with appro¬ 
priate and temporally adjacent master bias and flat field 
frames taken from the Gemini Science Archive. Each 
Gemini observation is the sum of six dithered 300 second 
integrations, stacked to eliminate cosmic rays. Gare was 
taken to avoid improperly rejecting pixels in the cores of 
bright stars when stacking the dithered frames. Though 
images A-G are well measured in each GMOS sub-image, 
the thick red-sensitive CCDs in use on GMOS-North are 
also excellent charged particle detectors, and the most ro¬ 
bust measurements thus come from stacked images with 
good cosmic ray rejection. 

The fluxes of SDSS J2222-I-2745 images A-C and ref¬ 
erence stars S1-S5 were all measured within a fixed aper¬ 
ture of diameter 1713, for the ALFOSC, MOSCA and 
GMOS data. The choice of aperture scale was motivated 
by containing most of the flux while avoiding contamina¬ 
tion from nearby sources. The effect on the time delay 
measurements of choosing an aperture a factor 2 smaller 
was tested, and found to produce changes in the time de¬ 
lay estimate within the uncertainties quoted below. The 
photometric uncertainty was determined by measuring 
the photometric scatter corresponding to the three in¬ 
dividual ALFOSC exposures and fitting the result as a 
function of the signal-to-noise of the quasar images. 

Figure [3] shows a histogram of seeing values measured 
in the individual ALFOSC monitoring exposures. In our 
further analysis, we excluded data points obtained when 
the FWHM seeing was worse than 1720, and we also 
excluded data where the sky background level was very 
high, e.g., caused by moonlight or twilight. This left 42 
distinct epochs of ALFOSC monitoring observations (not 
counting multiple exposures within the same night), with 
an average time span between observations of 14 days, 
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Fig. 1.— A single 600s ALFOSC exposure, with the locations of 
the five comparison stars S1-S5 indicated. 



Fig. 2. — Histogram of seeing values measured in individual AL¬ 
FOSC exposures. The dotted vertical line corresponds to the upper 
cutoff value for data used in our analysis. 

during the annual visibility period of this target. The 
photometric measurements used in our further analysis 
are listed in Table [H 

Significant flux variability of the three brightest images 
is evident even from a cursory inspection of the imaging 
data: Figure [3] qualitatively illustrates the flux changes 
over a 3.2 yr time span: In September 2011, the A and B 
images were both obviously brighter than the C image. 
By November 2014, the A and B components had both 
faded while approximately preserving the A:B flux ratio, 
while the C image had brightened to the point where it is 
clearly brighter than the B component and rivals the A 
component in magnitude. A quantitative representation 
of these trends is provided by the light curves of the three 
brightest components, shown in Figure 01 

3. THE A-B TIME DELAY 

Model predictions for the time delays between the six 
observed quasar images were reported in the work by 
Dahle et al. (2013). Being based on limited constraints 
from ground-based imaging (e.g., using only a single fam¬ 
ily of multiply lensed images, of the quasar itself), these 
predictions all have large uncertainty intervals of several 
hundred days. 
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Fig. 3. — Left Panel: MOSCA image from September 2011. Right 
Panel: ALFOSC image from November 2014, illustrating the flux 
variation of the C image relative to the A and B images of the 
quasar during our monitoring of SDSS J2222-1-2745. 
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Fig. 4.— Light curves of the quasar images A (blue symbols), B 
(green), and C (red). The light blue points in the upper part of the 
image represent the photometric measurements of the comparison 
stars S1-S5. 


The model prediction for the time delay between the 
A and B images reported by Dahle et al. (2013) is tab = 
112^158 days, where positive values imply that image 
A is leading image B. The time delay between A and 
C is predicted to be tqa = 1256^411 days. Given the 
much shorter predicted time delay between A and B, we 
would expect the quasar variability light curves of these 
two images to overlap over a large fraction of the time 
interval covered by our photometric monitoring. Hence, 
we focus initially on determining the value of tab- 
While a large number of different methods for deriv¬ 
ing the time delays of light curves of strong lenses have 
recently been proposed and tested (Liao et al. 2015), we 
focus here on two traditional methods for time delay mea¬ 
surements, detailed below. However, even prior to the 
implementation of these methods, an initial estimate of 
Tab of ~ 40 — 50 days was made by shifting and visually 
matching the light curves. Finding a plausible match was 
facilitated by the obvious overall variability trends seen 
in the A and B light curves, showing continuous declines 
over ^ 1.5 yr, followed by one observing season of nearly 
constant fluxes during 2014. Below, we use two different 
statistical methods to estimate tab in a more rigorous 
manner. 


3.1. minimization 

This method is based on shifting the two observed light 
curves A{ti) and cross-correlating each of them 
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Fig. 5. — The result of a minimization obtained by shifting and 
matching the observed light curves for A and B measured during 
the period September 2012 - August 2015. 

with a linearly interpolated light curve for the other im¬ 
age, and a(t), respectively), and finding the values 

of T and the magnitude offset Am which minimise the 
statistic 

X^(r) _ 1 [Mtj) - Ktj -t) + 

Adof(T) 2iVdof(r) ^ crid + o'M 

1 [a{ti +t) - B{U) + Am]'^ 

+ 2iVdof(T) 4" + ^ 


Fig. 6.— Dispersion calculated as a function of the shift tab 
between the two light curves. 

A ^-1 

E SkWkGk{Ck+i - Ckf 
D^{t) = min ^ - — -, (3) 

Am 

2 S SkWkGk 

k=l 

where Gk = 1 only if the measurements are from dif¬ 
ferent images and Gk =0 otherwise, Wk = 
are statistical weights, and 

Q _ / 1 l^fe+l ~ tk\ "G ^ (A\ 

\tk+i-tk\>5. 


Here, aA,i and aB,i are errors in the observed magni¬ 
tudes and aa,t and are errors in the interpolated mag¬ 
nitudes. We used data from all four seasons of system¬ 
atic monitoring and assumed that Auiab did not change 
between seasons (such changes might be produced by mi- 
crolensing). We did not interpolate between any set of 
points that spanned the gap between different seasons. 


The minimum value of 


X (r) 


as a function of the time 


Afdof(i-) 

delay is shown in Figure [5l We find a minimum value for 
a time delay tab = 42.01q 1 magnitude offset 

AmAB = 0-345tg;Qg5 (with uncertainties corresponding 


to Ax^ = 4). 


3.2. Dispersion method 

A shortcoming of the x^ minimisation method is the 
assumption of linear flux variations between measured 
points in the light curve. A more statistically robust al¬ 
ternative method for time delay estimates was proposed 
by Pelt et al. (1994; 1996). 

The method constructs a combined data set Gk(tk) 
from the two series Ai (with Na measurements) and 
Bj (with Nb measurements), sampling the quasar light 
curve: 


= { b'- Am Sl=g+r. (2) 

where fc = 1,..., A and N = Na + Nb ■ 

From this combined data set, the dispersion spectrum 
is defined by: 


We used 6 = 25 days as the decorrelation time scale, 
producing the dispersion spectrum shown in Figure |6l 
The minimum value is reached for tab = 47.7 ± 6.0 days 
and AmAB = 0.340 ± 0.007 (95% confidence intervals). 
The result has little sensitivity to the choice of <5; choos¬ 
ing a number twice or half the value quoted above 
changes the best fit values by a few percent, well within 
the uncertainty intervals quoted below. 

To estimate the uncertainties, we used a bootstrap pro¬ 
cedure. The light curve was smoothed using a running 
9-point median filter and the residuals of our data points 
were calculated with respect to the smoothed light curve. 
A set of 1000 simulated light curves were then gener¬ 
ated by randomly re-shuffling the residuals, and the er¬ 
rors were estimated from the spread of values of tab 
from the 1000 dispersion spectra. A combined A-l-B light 
curve G{tk) is represented by the union of the red and 
blue points plotted in Figure O In this figure, the ob¬ 
served light curve of the B component has been shifted by 
the Tab and AmAB values derived using the dispersion 
method. 


4. THE A-C TIME DELAY 

Having estimated tab , we next consider the time delay 
between the combined AB light curve and image C. For 
this purpose, we also include the 2009 photometry from 
SDSS DR8. All acceptable lens models of this system 
predict that image C leads all the other images, so we 
do not consider time delays t qa < 0. A x^ minimisa¬ 
tion using the procedure in § 13.11 yields a minimum at 
time delay tca = 717.4 ± 0.5 days, with uncertainties 
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Fig. 7.— The result of a minimization obtained by shifting 
and matching the observed light curve of C to the combined A+B 
light curve. 
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Fig. 8.— Dispersion calculated as a function of the shift tca 
between the observed light curve of image C and the combined 
A+B light curve. 

corresponding to Ax^ = 4. The magnitude offset is con¬ 
strained as ArriAC = 0.504 ± 0.013. The minimum value 

2 / N 

of + as a function of the time delay is shown in Fig¬ 
ure [T] No meaningful calculation could be done for time 
delay values greater than 830 days, since there would so 
far be minimal or no overlap between the monitoring 
light curves for such long time delays. 

We also esti mat ed tca using the dispersion method 
described in § 13.21 This produced the dispersion spec¬ 
trum shown in Figure [51 yielding a best-fit time delay 
of TCA = 722 ± 24 days and magnitude offset Atoac = 
0.483 ± 0.012 (95% confidence intervals). These values 
are consistent with the values produced by the min¬ 
imisation method. Since the dispersion method does not 
rely on assumptions about the quasar flux variation be¬ 
tween monitoring epochs, we consider the values and un¬ 
certainties produced by this method to be more reliable 
than those produced by interpolation and minimisa¬ 
tion. Hence, we adopt the time delay tca = 722 ±24 days 
in our further discussion. 

As shown in Figure I the best-fit value derived for 
TCA makes very distinctive predictions for the photo¬ 


metric behaviour of images A and B during the 2015 
and 2016 monitoring seasons: The A and B light curves 
are both predicted to reach a sharp minimum, spaced 
by their 47.7- day time delay, in mid-2015. This mini¬ 
mum corresponds to the lowest luminosity state of the 
quasar observed during the ~ 8 years of intrinsic quasar 
brightness variations spanned by all available photomet¬ 
ric data. Our most recent data points, from 2015 Aug. 
22, 24, and 26 UT (see Tabled]), show an increase in the 
flux of image A by ^ 0.1 magnitudes after passing the 
global minimum in the intrinsic quasar light curve, ex¬ 
actly as predicted by the best-fit tca value. Continued 
monitoring of images A and B, tracking their predicted 
further brightening during the fall of 2015, should sub¬ 
stantially reduce the size of the uncertainty intervals of 
both TAB and tca- 

The fluxes of A and B are predicted to further increase 
by > 0.5 magnitude during 2016. Observing these pre¬ 
dicted variations of image A and B through the end of 
2016 will strictly exclude any alternative time delay val¬ 
ues in the broad interval 830 < tca ^ 2100 days. We 
note that the upper limit of this interval extends well be¬ 
yond the range of time delays allowed by the lens model 
of Dahle et al. (2013). 

5. SUMMARY AND DISCUSSION 

Predictions from the lens model presented by Dahle 
et al. (2013) indicate that three independent time delays 
TAB) Ted and tfe in the system are all of order 100 days. 
The predicted time delays between image C and the pair 
(A-B) and between (A-B) and the central images (D-E- 
F) are of order 1300 days and 700-1000 days, respectively 
(with ~ 40% uncertainty). 

We have measured tab using two different methods 
and find results which are consistent at the 2a level. 
However, unlike the minimisation method, the dis¬ 
persion method does not make any assumptions about 
the photometric behaviour of the quasar images during 
the ^ 14 day sampling intervals, and we therefore adopt 
the dispersion method value of tab = 47.7 ± 6.0 days 
(95% confidence interval) as the most robust measure¬ 
ment. This is consistent with the model predicted value 
of TAB = 112+^^ days. The measured magnitude differ¬ 
ence Atuab = 0.340±0.007 is also consistent with the ra¬ 
tio of the model-predicted magnifications of /ta = 5-4lEg 

and hb = 4.7tE6 (Dahle et al. 2013). 

The measured tab time delay is very close to 1.5 syn¬ 
odic months, implying that it is possible to construct a 
continuous (apart from seasonal gaps) quasar light curve 
from the combination of data points from image A and 
image B which will be minimally affected by poor-S'/A 
measurements caused by bright moonlight. 

From our systematic monitoring campaign, there is so 
far only a year of overlap between the intrinsic quasar 
light curve derived from image C and those from A and B. 
However, the strong amplitude of the intrinsic brightness 
fluctuations of the quasar, the presence of earlier data 
points from 2009 and 2011, and the excellent match of 
the sharp flux minima observed in image C in 2013 and 
in images A and B in 2015 enable us to measure the tca 
time delay from the data set currently in hand. 

The best-fit value of tca = 722±24 days is in slight ten¬ 
sion with the model-predicted value of tca = 1256^411 
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days of Dahle et al. (2013). We note that ongoing work, 
refining the lens model by incorporating additional con¬ 
straints from newly identified multiply imaged source 
galaxies in this lens system, produce systematically lower 
model-predicted values for tca- The measured mag¬ 
nitude difference AmAC = 0.483 ±0.012 is consistent 
within the errors with the ratio of the model-predicted 
magnifications of and (Dahle 

et al. 2013). 

Our prediction of the timing of the expected 0.5 — 
0.6 magnitude rise in images A and B can be used to 
plan a more intensive campaign to monitor photomet¬ 
ric and spectroscopic changes in the quasar during this 
phase in 2016. This would e.g., allow a reverberation 
mapping study of the z = 2.82 quasar. 

A subsequent corresponding brightness increase should 
occur for the three fainter images around 2018, given the 
Dahle et al. (2013) prediction for the time delays be¬ 
tween the (A-B) pair and the (D-F) trio. During 2015, 
these fainter images should be entering the decline phase 
seen in images A and B in 2012-2013, before reaching the 
plateau seen in 2014. Provided that even rough estimates 
of the time delays can be made from observations in the 
decline and subsequent plateau phase, a targeted cam¬ 
paign at the predicted time of the subsequent rise could 
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provide accurate time delay estimates for the three de- 
magnified images in the cluster center. 

We have recently obtained 5-band HST optical/NIR 
imaging (GO-13337; PI: Sharon) of SDSS J2222±2745 
which confirms the existence of the sixth quasar image 
and reveals additional multiply lensed background galax¬ 
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model of this system and provide updated predictions 
for the unmeasured time delays. The HST data will also 
be crucial for modelling the brightness distribution of 
the brightest galaxies in the cluster core and removing 
their contribution to the measured flux of images D-F in 
ground-based photometric monitoring data. 

To summarise: The steep quasar brightness fluctua¬ 
tions evident from the light curve of the leading C image 
provide strong constraining power to measurements of 
the time delays between all six known images in this lens 
system. The range of clustercentric radii over which such 
measurements are possible is truly unique for this system. 
The next 2-3 observing seasons will be instrumental in 
pinning down these time delays, enabling a wide range of 
physical studies of both the cluster lens and the quasar 
source. 
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Fig. 9.— The combined A (open blue circles) + B (open green circles) light curve, together with the C (filled red circles) light curve, the 
latter shifted by = 722 days and AmAC = 0.483. The dates and magnitude values correspond to the light curve of the A image. 
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TABLE 1 

g-BAND PHOTOMETRY OF THE THREE BRIGHTEST QUASAR IMAGES IN SDSS 

J2222+2745 


MJD 

Image A 

Image B 

Image C 

Observatory 

Detector 

5122.200 

21.100± 

0.030 

21.430± 

0.040 

21.510± 

0.050 

SDSS 

SDSS 

5828.935 

20.949± 

0.009 

21.422± 

0.009 

21.987± 

0.009 

NOT 

MOSCA 

6182.963 

21.092± 

0.009 

21.441± 

0.009 

22.111± 

0.009 

NOT 

MOSCA 

6185.985 

21.094± 

0.009 

21.465± 

0.009 

22.111± 

0.009 

NOT 

MOSCA 

6244.907 

21.170± 

0.014 

21.473± 

0.017 

22.145± 

0.030 

NOT 

ALFOSC 

6244.911 

21.148± 

0.014 

21.524± 

0.014 

22.083± 

0.019 

NOT 

ALFOSC 

6244.915 

21.155± 

0.014 

21.507± 

0.014 

22.100± 

0.021 

NOT 

ALFOSC 

6266.907 

21.187± 

0.014 

21.509± 

0.015 

22.101± 

0.023 

NOT 

ALFOSC 

6266.911 

21.214± 

0.014 

21.513± 

0.015 

22.080± 

0.022 

NOT 

ALFOSC 

6270.830 

21.218± 

0.015 

21.529± 

0.014 

22.122± 

0.021 

NOT 

ALFOSC 

6270.834 

21.213± 

0.015 

21.518± 

0.014 

22.089± 

0.020 

NOT 

ALFOSC 

6270.838 

21.224± 

0.015 

21.499± 

0.014 

22.070± 

0.020 

NOT 

ALFOSC 

6291.811 

21.301± 

0.019 

21.527± 

0.024 

22.120± 

0.037 

NOT 

ALFOSC 

6291.815 

21.242± 

0.015 

21.516± 

0.017 

22.168± 

0.028 

NOT 

ALFOSC 

6291.819 

21.288± 

0.014 

21.568± 

0.015 

22.201± 

0.024 

NOT 

ALFOSC 

6430.182 

21.396± 

0.015 

21.792± 

0.015 

22.189± 

0.019 

NOT 

ALFOSC 

6430.189 

21.441± 

0.015 

21.815± 

0.015 

22.243± 

0.019 

NOT 

ALFOSC 

6430.197 

21.413± 

0.015 

21.741± 

0.014 

22.198± 

0.019 

NOT 

ALFOSC 

6446.168 

21.363± 

0.015 

21.683± 

0.020 

22.191± 

0.028 

NOT 

ALFOSC 

6446.182 

21.327± 

0.014 

21.655± 

0.019 

22.111± 

0.025 

NOT 

ALFOSC 

6454.178 

21.430± 

0.014 

21.821± 

0.019 

22.234± 

0.026 

NOT 

ALFOSC 

6454.182 

21.443± 

0.015 

21.763± 

0.018 

22.261± 

0.026 

NOT 

ALFOSC 

6454.186 

21.382± 

0.014 

21.764± 

0.018 

22.233± 

0.025 

NOT 

ALFOSC 

6457.183 

21.388± 

0.007 

21.767± 

0.014 

22.254± 

0.015 

NOT 

ALFOSC 

6457.190 

21.401± 

0.011 

21.800± 

0.014 

22.260± 

0.016 

NOT 

ALFOSC 

6457.198 

21.373± 

0.011 

21.778± 

0.014 

22.260± 

0.016 

NOT 

ALFOSC 

6477.090 

21.452± 

0.015 

21.765± 

0.014 

22.207± 

0.019 

NOT 

ALFOSC 

6477.097 

21.455± 

0.014 

21.749± 

0.014 

22.251± 

0.018 

NOT 

ALFOSC 

6477.105 

21.463± 

0.014 

21.767± 

0.014 

22.295± 

0.018 

NOT 

ALFOSC 

6488.130 

21.444± 

0.010 

21.721± 

0.015 

22.220± 

0.016 

NOT 

ALFOSC 

6488.137 

21.472± 

0.010 

21.728± 

0.015 

22.249± 

0.016 

NOT 

ALFOSC 

6488.145 

21.443± 

0.010 

21.719± 

0.015 

22.227± 

0.016 

NOT 

ALFOSC 

6507.150 

21.470± 

0.016 

21.781± 

0.021 

22.289± 

0.030 

NOT 

ALFOSC 

6507.158 

21.455± 

0.016 

21.757± 

0.021 

22.265± 

0.030 

NOT 

ALFOSC 

6507.165 

21.445± 

0.015 

21.751± 

0.020 

22.283± 

0.030 

NOT 

ALFOSC 

6518.141 

21.451± 

0.015 

21.761± 

0.015 

22.288± 

0.020 

NOT 

ALFOSC 

6518.148 

21.471± 

0.014 

21.808± 

0.014 

22.302± 

0.018 

NOT 

ALFOSC 

6518.156 

21.459± 

0.013 

21.789± 

0.014 

22.281± 

0.017 

NOT 

ALFOSC 

6533.039 

21.475± 

0.015 

21.812± 

0.014 

22.213± 

0.017 

NOT 

ALFOSC 

6533.050 

21.502± 

0.015 

21.824± 

0.015 

22.262± 

0.019 

NOT 

ALFOSC 

6533.058 

21.507± 

0.015 

21.820± 

0.015 

22.235± 

0.019 

NOT 

ALFOSC 

6533.065 

21.466± 

0.014 

21.837± 

0.015 

22.285± 

0.022 

NOT 

ALFOSC 

6543.103 

21.497± 

0.014 

21.799± 

0.015 

22.212± 

0.020 

NOT 

ALFOSC 

6543.110 

21.496± 

0.014 

21.830± 

0.015 

22.208± 

0.019 

NOT 

ALFOSC 

6543.118 

21.501± 

0.014 

21.788± 

0.015 

22.187± 

0.019 

NOT 

ALFOSC 

6549.047 

21.528± 

0.015 

21.836± 

0.020 

22.138± 

0.026 

NOT 

ALFOSC 

6549.054 

21.461± 

0.014 

21.815± 

0.017 

22.220± 

0.023 

NOT 

ALFOSC 

6557.850 

21.555± 

0.014 

21.811± 

0.017 

22.197± 

0.023 

NOT 

ALFOSC 

6557.854 

21.552± 

0.014 

21.853± 

0.017 

22.137± 

0.022 

NOT 

ALFOSC 

6557.858 

21.523± 

0.015 

21.826± 

0.018 

22.172± 

0.024 

NOT 

ALFOSC 

6565.970 

21.573± 

0.012 

21.800± 

0.015 

22.169± 

0.014 

NOT 

ALFOSC 

6565.977 

21.559± 

0.013 

21.829± 

0.014 

22.196± 

0.015 

NOT 

ALFOSC 

6565.985 

21.531± 

0.011 

21.811± 

0.015 

22.137± 

0.014 

NOT 

ALFOSC 

6573.944 

21.579± 

0.010 

21.810± 

0.015 

22.162± 

0.014 

NOT 

ALFOSC 

6573.951 

21.560± 

0.010 

21.806± 

0.015 

22.182± 

0.014 

NOT 

ALFOSC 

6573.958 

21.579± 

0.009 

21.853± 

0.015 

22.170± 

0.014 

NOT 

ALFOSC 

6578.051 

21.520± 

0.018 

21.737± 

0.023 

22.088± 

0.031 

NOT 

ALFOSC 

6578.055 

21.542± 

0.020 

21.806± 

0.025 

22.078± 

0.031 

NOT 

ALFOSC 

6578.059 

21.577± 

0.017 

21.841± 

0.022 

22.149± 

0.027 

NOT 

ALFOSC 

6602.891 

21.634± 

0.015 

21.868± 

0.014 

22.135± 

0.015 

NOT 

ALFOSC 

6602.899 

21.621± 

0.014 

21.889± 

0.015 

22.160± 

0.018 

NOT 

ALFOSC 

6602.906 

21.585± 

0.011 

21.862± 

0.015 

22.156± 

0.014 

NOT 

ALFOSC 

6618.884 

21.650± 

0.014 

21.899± 

0.015 

22.179± 

0.018 

NOT 

ALFOSC 

6618.888 

21.638± 

0.014 

21.880± 

0.015 

22.163± 

0.017 

NOT 

ALFOSC 

6618.891 

21.602± 

0.014 

21.903± 

0.016 

22.139± 

0.020 

NOT 

ALFOSC 

6788.186 

21.572± 

0.015 

21.916± 

0.019 

22.160± 

0.022 

NOT 

ALFOSC 

6788.193 

21.597± 

0.015 

22.003± 

0.020 

22.133± 

0.022 

NOT 

ALFOSC 

6815.173 

21.556± 

0.015 

21.942± 

0.015 

22.113± 

0.016 

NOT 

ALFOSC 

6815.181 

21.570± 

0.015 

21.960± 

0.015 

22.112± 

0.016 

NOT 

ALFOSC 

6815.188 

21.565± 

0.015 

21.942± 

0.015 

22.101± 

0.016 

NOT 

ALFOSC 

6835.176 

21.614± 

0.014 

21.927± 

0.014 

22.076± 

0.015 

NOT 

ALFOSC 

6835.184 

21.584± 

0.015 

21.914± 

0.014 

22.080± 

0.015 

NOT 

ALFOSC 

6835.191 

21.571± 

0.015 

21.924± 

0.015 

22.117± 

0.016 

NOT 

ALFOSC 

6847.162 

21.584± 

0.014 

21.901± 

0.015 

22.022± 

0.016 

NOT 

ALFOSC 

6847.169 

21.574± 

0.014 

21.916± 

0.015 

22.062± 

0.016 

NOT 

ALFOSC 
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TABLE 1 — Continued 


MJD 

Image A 

Image B 

Image C 

Observatory 

Detector 

6847.179 

21.589± 

0.014 

21.908± 

0.016 

22.060± 

0.016 

NOT 

ALFOSC 

6864.105 

21.584± 

0.015 

21.959± 

0.015 

22.017± 

0.015 

NOT 

ALFOSC 

6874.056 

21.610± 

0.014 

21.914± 

0.015 

21.950± 

0.015 

NOT 

ALFOSC 

6874.063 

21.612± 

0.015 

21.934± 

0.014 

21.952± 

0.014 

NOT 

ALFOSC 

6894.020 

21.610± 

0.012 

21.937± 

0.014 

21.963± 

0.014 

NOT 

ALFOSC 

6894.028 

21.592± 

0.014 

21.902± 

0.014 

21.903± 

0.014 

NOT 

ALFOSC 

6894.035 

21.599± 

0.015 

21.892± 

0.014 

21.943± 

0.014 

NOT 

ALFOSC 

6918.096 

21.565± 

0.014 

21.956± 

0.018 

21.915± 

0.017 

NOT 

ALFOSC 

6928.011 

21.619± 

0.015 

21.952± 

0.014 

21.956± 

0.014 

NOT 

ALFOSC 

6928.019 

21.620± 

0.014 

21.987± 

0.015 

21.944± 

0.015 

NOT 

ALFOSC 

6928.026 

21.625± 

0.013 

21.942± 

0.014 

21.959± 

0.014 

NOT 

ALFOSC 

6944.979 

21.600± 

0.014 

21.933± 

0.015 

21.923± 

0.014 

NOT 

ALFOSC 

6944.986 

21.628± 

0.015 

21.921± 

0.014 

21.960± 

0.014 

NOT 

ALFOSC 

6944.993 

21.625± 

0.014 

21.933± 

0.015 

21.900± 

0.015 

NOT 

ALFOSC 

6973.943 

21.629± 

0.018 

21.931± 

0.024 

21.699± 

0.019 

NOT 

ALFOSC 

6973.950 

21.627± 

0.018 

21.920± 

0.024 

21.726± 

0.019 

NOT 

ALFOSC 

6973.958 

21.524± 

0.019 

21.952± 

0.027 

21.709± 

0.021 

NOT 

ALFOSC 

6986.896 

21.643± 

0.015 

21.883± 

0.018 

21.654± 

0.015 

NOT 

ALFOSC 

6986.903 

21.642± 

0.014 

21.920± 

0.017 

21.683± 

0.014 

NOT 

ALFOSC 

6986.911 

21.623± 

0.014 

21.918± 

0.017 

21.668± 

0.014 

NOT 

ALFOSC 

7141.226 

21.651± 

0.031 

22.003± 

0.039 

21.788± 

0.034 

NOT 

ALFOSC 

7143.207 

21.696± 

0.024 

22.050± 

0.030 

21.772± 

0.024 

NOT 

ALFOSC 

7143.211 

21.744± 

0.020 

22.067± 

0.027 

21.809± 

0.022 

NOT 

ALFOSC 

7143.214 

21.696± 

0.018 

21.966± 

0.026 

21.742± 

0.019 

NOT 

ALFOSC 

7168.539 

21.789± 

0.009 

22.031± 

0.009 

21.836± 

0.009 

Gemini 

CMOS 

7169.188 

21.750± 

0.014 

21.986± 

0.016 

21.822± 

0.014 

NOT 

ALFOSC 

7169.195 

21.775± 

0.014 

22.061± 

0.016 

21.824± 

0.014 

NOT 

ALFOSC 

7169.203 

21.800± 

0.014 

22.014± 

0.017 

21.826± 

0.014 

NOT 

ALFOSC 

7169.518 

21.787± 

0.009 

22.025± 

0.009 

21.830± 

0.009 

Gemini 

CMOS 

7191.154 

21.771± 

0.014 

22.044± 

0.015 

21.818± 

0.014 

NOT 

ALFOSC 

7191.161 

21.777± 

0.014 

22.029± 

0.016 

21.822± 

0.014 

NOT 

ALFOSC 

7191.168 

21.783± 

0.014 

21.996± 

0.016 

21.857± 

0.014 

NOT 

ALFOSC 

7196.183 

21.784± 

0.015 

22.052± 

0.014 

21.828± 

0.014 

NOT 

ALFOSC 

7196.190 

21.809± 

0.014 

22.039± 

0.015 

21.817± 

0.014 

NOT 

ALFOSC 

7196.197 

21.807± 

0.014 

22.059± 

0.015 

21.819± 

0.014 

NOT 

ALFOSC 

7218.080 

21.794± 

0.016 

22.110± 

0.022 

21.740± 

0.015 

NOT 

ALFOSC 

7218.087 

21.763± 

0.016 

22.054± 

0.021 

21.731± 

0.016 

NOT 

ALFOSC 

7228.476 

21.791± 

0.009 

22.096± 

0.009 

21.793± 

0.009 

Gemini 

CMOS 

7256.088 

21.687± 

0.014 

22.066± 

0.016 

21.748± 

0.014 

NOT 

ALFOSC 

7258.077 

21.725± 

0.014 

22.079± 

0.014 

21.746± 

0.014 

NOT 

ALFOSC 

7258.084 

21.697± 

0.015 

22.121± 

0.015 

21.723± 

0.015 

NOT 

ALFOSC 

7258.091 

21.696± 

0.015 

22.066± 

0.015 

21.733± 

0.015 

NOT 

ALFOSC 

7260.148 

21.691± 

0.014 

22.087± 

0.017 

21.763± 

0.014 

NOT 

ALFOSC 


Note. — The Modified Julian Days (MJD) column gives the date of the observation relative to MJD= 50000. 




